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ABSTRACT 

We present i7-band VLT/NACO polarized light images of the Herbig Ae/Be star HD169142 probing 
its protoplanetary disk as close as ^0.1" to the star. Our images trace the face-on disk out to ~1.7" 
(^-^250 AU) and reveal distinct sub-structures for the first time: 1) the inner disk (<20 AU) appears 
to be depleted in scattering dust grains; 2) an unresolved disk rim is imaged at ^25 AU; 3) an annular 
gap extends from ~40 - 70 AU; 4) local brightness asymmetries are found on opposite sides of the 
annular gap. We discuss different explanations for the observed morphology among which ongoing 
planet formation is a tempting - but yet to be proven - one. Outside of ^85 AU the surface brightness 
drops off roughly cx r~^^^, but describing the disk regions between 85-120 AU / 120-250 AU separately 
with power-laws cx r~^'^/cx r~^'^ provides a better fit hinting towards another discontinuity in the disk 
surface. The flux ratio between the disk integrated polarized light and the central star is ~ 4.1 • 10^"^. 
Finally, combining our results with those from the literature, ~40% of the scattered light in the H- 
band appears to be polarized. Our results emphasize that IID169142 is an interesting system for 
future planet formation or disk evolution studies. 

Subject headings: stars: pre-main sequence — stars: formation — protoplanetary disks — planet-disk 
interactions — stars: individual (HD169142) 



1. INTRODUCTION 

To study the physical and chemical conditions for 
planet formation and to search for morphological evi- 
dence of forming planets, the inner few tens of AU of pro- 
toplanetary disks have to be investigated. In particular, 
annular gaps in those disks are considered to be possi- 
ble signatures of ongoing planet formation. In this letter 
we report the detection of such gaps in the disk around 
the Herbig Ae/Be star HD169142 using polarimetric dif- 
ferential imaging (PDI) . PDI is a powerful high-contrast 
technique to study the dusty surface layer of protoplane- 
tary disks very close to the star (e.g., Quanz et al. 2011; 
Hashimoto et al. 2011; Muto et al. 2012; Quanz et al. 
2012; Hashimoto et al. 2012; Tami et al. 2012; Kusakabe 
et al. 2012; Mayama et al. 2012; Grady et al. 2013). 

Our new i7-band PDI observations of HD169142 re- 
solve for the first time distinct structures in the inner 
disk regions with high signal-to-noise and allow an in- 
depth analysis of the disk profile and morphology. The 
basic stellar parameters of HD169142 are summarized in 
Table 1. 

HD169142 was classified as a Herbig group lb member 
(Meeus et al. 2001), indicating that its spectral energy 
distribution (SED) rises in the mid-infrared (MIR), but 
does not show a 10 /im silicate emission feature (e.g., 
van Boekel et al. 2005). However, PAH emission fea- 
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TABLE 1 
Basic parameters of HD169142. 



Parameter 


Value for HD169142 


Reference'' 


RA (J2000) 


18''24™■29^79 


(1) 


DEC (J2000) 


-29°46'49".22 


(1) 


J [mag] 


7.31±0.02 


(1) 


H [mag] 


6.91±0.04 


(1) 


Ks [mag] 


6.41±0.02 


(1) 


Sp. Type 


A9III/IVe / A7V 


(2),(3) 


Dsini [kms"-'] 


55±5 


(2) 


Age [Myr] 


1-5 / 12 / 3-12 


(2),(3),(4) 


[Fe/H] 


-0.5 ±0.1 / -0.25 - -0.5 


(2),(5) 


log 9 


3.7±0.1 / 4.0-4.1 


(2),(5) 


T,a [K] 


7500±200 / 6500 / 7650±150 


(2),(3),(5) 


Mass [Mq] 


~1.65 


(3) 


R, [Rq] 


~1.6 


{3),(5) 


L* [Lq] 


~8.6 


(3) 


M [10-9 M0yr-l] 


~3.1 / <1.25±0.55 


(3),(4) 


Distance [pc] 


145*^ / 151 


(6),(3) 



''References — (1) 2MASS point source catalog (Cutri et al. 
2003), (2) Guimaraes et al. (2006), (3) Blondel & Djie (2006), 
(4) Grady et al. (2007), (5) Meeus et al. (2010), (6) Sylvester et al. 
(1996). 

''adopted value in this paper 



tures have been detected in MIR spectroscopy (Meeus 
et al. 2001; van Boekel et al. 2005) and were spatially 
resolved at 3.3 /im (Habart et al. 2006). The disk was 
previously resolved in polarized light in the NIR from the 
ground (Kuhn et al. 2001; Hales et al. 2006), but those 
data did not allow an in-depth analysis of the disk struc- 
ture and profile. Grady et al. (2007) and Fukagawa et al. 
(2010) detected the disk in scattered light at 1.1 //m with 
HST/NICMOS coronagraphy and in H-hand from the 
ground, respectively, and they derived azimuthally av- 
eraged surface brightness profiles. In all direct imaging 
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TABLE 2 

Summary of NACO/PDI observations of HD169142 on July 25, 2012. 



Filter 


DIT 


X NDIT-'^ 


Dither 


Average 


Average 


Average (to)'* 








l)(>si(i(>lis'' 


aiiiiiass 


StH'illg'' 




iVBl.64 


1 


s X 15 


3 


~1.02 


~ 1.57" 


~ 14 ms 


H 


1 


s X 45 


2x6 


~1.06 


~ 1.04" 


~ 22 ms 


NBl.M 


1 


s X 20 


3 


~1.11 


~ 1.07" 


~ 20 ms 



''Detector integration time (DIT) X number of integrations (NDIT), i.e., total 
integration time per dither position and per retarder plate position. 

''At each dither position 4 diflferent HWP positions (0.0°, -22.5°, -45.0°, - 
67.5°) were used. 

"^DIMM seeing in the optical as measured from seeing monitor . 
Average coherence time of the atmosphere calculated by the Real Time 
Computer of the AO system. 



studies the disk was eonsistent with being seen faee-on 
and no sub-structures were detected. In the MIR the 
disk was resolved at 12, 18, 18.8 and 24.5 /zm (Honda 
et al. 2012; Marinas et al. 2011) and most of this emis- 
sion is thought to arise from disk regions between ~30 - 
60 AU. From dust continuum and molecular line observa- 
tions the mass and outer radius of the disk are estimated 
to be -0.005 - 0.04 Mq and ~235 AU (Raman et al. 
2006; Dent et al. 2006; Panic et al. 2008; Meeus et al. 
2010; Sandell et al. 2011). The kinematic pattern and 
molecular line profiles are best fitted with a disk incli- 
nation of ~13° (Raman et al. 2006; Panic et al. 2008). 
The lack of a silicate emission feature at 10 /ini can be 
explained with a lack of small (<3-5 /um) dust grains in 
the inner ~ 10 - 20 AU of the disk possibly due to grain 
growth (van Boekel et al. 2005). From SED modeling a 
disk hole was inferred for the inner disk regions (Grady 
et al. 2007; Meeus et al. 2010; Honda et al. 2012). 

2. OBSERVATIONS AND DATA REDUCTION 

The observations were carried out with VLT/NACO 
(Lenzen et al. 2003; Rousset et al. 2003) in the H and 
NBl.QA filter. A description of the camera properties 
and the PDI observing mode is given in Quanz et al. 
(2011). In short, a WoUaston prism splits the light in an 
ordinary and extraordinary beam with orthogonal, linear 
polarization directions. Both beams are imaged simul- 
taneously on the detector. The polarization direction 
can be changed by rotating a half- wave plate (HWP). 
One full polarization cycle includes exposures at four 
HWP positions (0° and -45° for Stokes Q, -22.5° and 
—67.5° for Stokes U). We summarize the observations 
for HD169142 in Table 2. The central few pixels of the 
point spread function (PSF) were saturated in the H fil- 
ter. The exposures in the NBl.&A filter were unsaturated 
and used for photometric calibration (see below). 

NACO suffers from instrumental polarization and 
cross-talk between the different Stokes components 
(Witzel et al. 2011; Quanz et al. 2011). To obtain re- 
liable results in PDI observations these effects need to 
be corrected for. Having obtained data sets for five more 
Herbig Ae/Be stars allowed us to do a systematic anal- 
ysis and we refer to an upcoming ptiblication for a de- 
tailed description of our data reduction and calibration 
approach (Avenhaus et al., in prep.). The key steps are: 

All exposures were dark current corrected and flat- 
fielded. Bad pixels were replaced with the mean value 
of surrounding pixels. Each quadrant of the detector 



showed a noise pattern affecting every other row, which 
was eliminated by subtracting the mean value of each row 
from each pixel in that row. Detector regions close to the 
star were excluded to comptite the mean values. The or- 
dinary and extraordinary images in each exposure were 
extracted and the stellar position in all individual images 
determined. We rebinned the images to a factor three 
higher resolution (bicubic interpolation) and shifted all 
images to a common center (bilinear interpolation). As- 
suming that the central star is unpolarized in if-band 
we corrected instrumental polarization and polarization 
cross-talk effects^. For each dither position the fractional 
Stokes parameters pq and pu were computed using the 
double ratio approach (see e.g., Tinbergen 1996; Schmid 
et al. 2006): 



PQ 



Rq + 1 



Pu 



Rc/- 1 
Rc/ + 1 



(1) 



with 



Rq 



Old/ extra 
r- 45° IT-i5° 
ord / extra 



Rl 



-22.5° 
rd 



22.5 
extra 



7—67.5° 
-'ord 



/ 7—67.5° 
/ -'extra 



(2) 

Here, the subscripts refer to either the ordinary or ex- 
traordinary beam and the superscripts refer to the angu- 
lar position of the HWP. The final pq and pu images were 
obtained by averaging over all dither positions. From 
PQ and Pu and the intensity images / we obtained the 
Stokes Q and U images, from which we computed the 
radial Stokes parameters and U^: 



Qr = +Q cos 2(j)+U sin 2(/) 
= -Q sm2(j) + U cos 2(/) 



with 



arctan 



y-yo 



(3) 
(4) 

(5) 



being the polar angle of a given position (x, y) on the 
detector, (xq, i/o) denoting the central position of the star 
(see, e.g., Schmid et al. 2006) and 9 being an offset due to 
polarization cross-talk effects (Avenhaus et al., in prep.), 
which amounts to 7° in our case. Qr is equivalent to the 
polarized flux P under the assumption that the polarized 
flux has only a tangential component, but, in contrast to 

^ The results were consistent with those using the approach de- 
scribed in Quanz et al. (2011). 



Gaps in the HD169142 protoplanetary disks revealed by polarimetric imaging 



3 



-150 



0.5" 



-0.5" 



Distance (AU) 
-100 -50 50 100 150 



^^^^ 



Distance (arcsec) 
0.5" 0" 0.5" 



■150 150 



100 100 



50 S 
< 



g 

s 

-50 b -50 -K 

-100 -100-,* 

-150 -150 




0.5" 



-0.5' 



-0.5" 0" 0.5" 1" 
Distance (arcsec) 



-150 -100 -50 50 100 150 
Distance (AU) 




-100 -50 50 100 150 
Distance (AU) 



Degrees east of north 



Fig. 1. — NACO/PDI observations of HD169142 in the H band, a) Final Qr image scaled with to compensate for the decrease in 
stellar flux (image shown in a linear stretch). The position of the central star is indicated by the red cross. Saturated pixels in the central 
regions have been masked out. Our data reveal a bright inner ring, a large gap and a smooth outer disk in polarized light. A brightness 
dip in the ring and a residual AO feature are indicated by arrows, b) Ur with the same scaling and stretch as for the Qr image, c) Intensity 
image also scaled with r^. Features from the AO system and the telescope spiders are clearly seen, d) Polar coordinate mapping of Qr- 
The innermost, masked out region is less than 0.1" in diameter. The red line traces the peak brightness of the inner ring. 



P, Qr is free of any biases introduced from computing 
the squares of the Q and U components (Schmid et al. 
2006). Under the same assumption, and assuming the 
correction for instrumental effects was done properly, [/j. 
only contains noise and can be used to estimate the error 
in Qr via 



AQ, 



(6) 



Here, AQ^ refers to the la uncertainty in Qr, crfj is the 
variance in the final image and n^os the number of 
resolution elements in the region of interest. 

The calibration of the surface brightness level in the 
final Qi- image was done using the non-saturated images 
in the NBl.GA filter following the description in Quanz 
et al. (2011). We estimate that the absolute fiux calibra- 
tion is good to ~30%. 

3. RESULTS AND ANALYSIS 



Figure 1 shows the final Qr and images, an intensity 
image, and a polar coordinate mapping of the Q,- image. 
While the Qr image shows extended polarized flux that 
we interpret as the disk around HD169142, the im- 
age does not reveal any significant signal as expected for 
scattering of dust particles without any preferred align- 
ment. 

Figure 1 a) reveals a protoplanetary disk with a com- 
plex radial morphology. Very close to the star, but out- 
side the innermost, saturated pixels, the detected polar- 
ized flux is low and it increases going outwards to reach 
a ring-like maximum at ~0.17" (~25 AU). Outside of 
the narrow ring the surface brightness decreases and an 
annular gap stretches from ~0.28-0.48" (--40 - 70 AU) 
with a local minimum at ~0.38" (~55 AU). Outside of 
~0.52" (~75 AU) the disk surface brightness drops off 
smoothly. 

The bright inner ring does not appear fully circular 
(Figure 1 panels a and d) and the radius of the peak 
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brightness varies between ~0.16" and ~0.19" with a 
mean value of ^0.17". A dip in the brightness of the 
ring seems to be present at position angle (PA) «80° 
(east of north). On the opposite side of the annular gap 
at PA«90°, next to an AO artifact, a region of enhanced 
brightness is seen stretching a bit into the gap. The full- 
width-half-maximum of the ring varies slightly with PA 
but the ring is not spatially resolved. Using the ring we 
estimated an upper limit in disk inclination of i < 20° , so 
that i is still in agreement with the fluctuations in ring 
radius. However, our image is consistent with a face-on 
disk. Analyzing the mean surface brightness as a func- 
tion of azimuth for different radii also supports a small 
inclination angle as no significant variation is seen (Muto 
et al. (2012) used this method for the SAO206462 disk). 

In Figure 2 we show the azimuthally averaged surface 
brightness profile and the signal-to-noise profile for the 
Qr and images. Fitting a power-law to the bright- 
ness profile between 85-250 AU yields an exponent of 
— 3.31±0.11, where the error denotes the 95% confidence 
level (black dashed line in top panel of Figure 2). How- 
ever, fitting the regions between 85-120 AU and 120-250 
AU separately with power-laws yields a mathematically 
better fit in a sense (black solid line). In this case the 
power-law exponents are — 2.641q and — 3.90j^Q |^° for 
the inner and outer region, respectively. The power-law 
fit to the inner region connects well with the bright inner 
ring at ~25 AU. A 0.15"x0.1" region centered on the 
AO feature seen in Figure 1 has been masked out before 
the power-law fits were computed. 

The maximum surface brightness in polarized fiux of 
S ~ 11.3 mag/arcsec^ is reached at the location of the 
bright ring. Close to the inner edge of the outer disk 
(-75 AU) we find ^14.3 mag/arcsec^ and at -220 AU 
wl8.1 mag/arcsec^. Inside the gap (—55 AU) the sur- 
face brightness is reduced by a factor of —3 compared to 
the corresponding value of the power- law fit. Integrating 
the polarized surface brightness between —14.5-250 AU 
and comparing it to the iJ-band fiux of the central star 
yields a flux ratio of F^^^^"""""^ / Fst^r ~ 4.1 • lO'^. Here, 
the regions between —14.5-40 AU (the ring), —40-70 AU 
(the gap), and —70-250 AU (the outer disk) contribute 
-49%, -14%, and -37% to the integrated disk flux, re- 
spectively. 

4. DISCUSSION 

The face-on orientation and general extent of the disk 
have been reported in earlier scattered or polarized light 
studies (Kuhn et al. 2001; Hales et al. 2006; Grady et al. 
2007; Fukagawa et al. 2010). However, no sub-structures 
were revealed so far. 

4.1. The bright inner ring 

We interpret the bright ring as the inner rim of the 
main disk. Its location flts well to recent model predic- 
tions (Meeus et al. 2010; Honda et al. 2012) but is not in 
agreement with a rim at 44 AU as suggested by Grady 
et al. (2007). The drop in polarized light inside of the rim 
is significant and supports the idea that a hole might be 
present in the inner <20 AU. The brightness dip in the 
rim at PA«80° hints to additional substructure or local 
change in dust grain properties. Although local minima 
in polarized flux can result from the scattering function 
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Fig. 2. — Azimuthally averaged surface brightness profile and 
SNR for the HD169142 disk observed with NACO/PDI in the H 
band. Top: Surface brightness in polarized flux. Error bars are Icr 
and include the the error defined in eq. (6) and the standard devi- 
ation from the mean flux in each annulus. The triangle represents 
a 1(7 upper limit. In the gray shaded region the noise dominates. 
The blue dash-dotted curve shows Scr the detection limit. The 
black dashed line is a power-law fit to the brightness profile be- 
tween 85-250 AU. The black solid line is a broken power-law with 
a break point at 120 AU (see text). Middle: The observed sur- 
face brightness profile (red data points) and both power-law fits 
(dashed line = single power-law; broken power-law = solid line) 
(divided by the single power-law fit. A broken power-law clearly 
provides a better fit to the data. Bottom: SNR for the Qr (red 
stars) and Ur (blue crosses) images. The signal in Ur is consistent 
with random noise, while in Qr the disk can clearly be traced out 
to ~1.7". 

of the dust grains even for small disk inclinations (see, 
e.g., Perrin et al. 2009; Muto et al. 2012) also the outer 
disk regions should show a minimum at the same PA, 
which is not the case. 

4.2. The annular gap 

To our knowledge this is the first time that such a sym- 
metric and finite gap-like structure has been detected on 
the surface layer of a protoplanetary disk. This gap does 
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not appear to be completely void of scattering material as 
the remaining surface brightness is above our detection 
limits. Seeing a gap-like feature in scattered light could 
indicate a deficiency in surface density possibly penetrat- 
ing deeper into the disk mid-plane. If so, one (or several) 
forming planet (s) could be responsible for the observed 
morphology. Theory and simulations predict that plan- 
ets open annular gaps in disks and the gap width and 
depth depend on the planet mass and on local disk prop- 
erties (e.g., Lin & Papaloizou 1986; Bryden ct al. 1999; 
Ruge et al. 2013). Jang-Condell & Turner (2012) simu- 
lated how planet-induced gaps would alter scattered light 
images of disks. They showed that depending on the 
mass of the planet the gap would not necessarily appear 
empty (cf. Dong et al. 2012). In addition, SPH simula- 
tions suggest that localized streams through the disk gap 
connect the outer and inner disk with the forming planet 
(e.g., Ruge et al. 2013). One could speculate whether 
the observed local brightness asymmetries on opposite 
sides of the gap (see section 3) could be signposts of such 
features. Finally, according to the simulations by Jang- 
Condell & Turner (2012), the impact of gaps and plan- 
ets on the SEDs of star+disk systems is limited. This 
could explain why previous SED modeling attempts for 
HD169142 did not reveal this feature. 

One could in principle infer mass constraints for the 
companion from the geometry of the gap (Lin & Pa- 
paloizou 1986; Bryden et al. 1999). However, this ex- 
ercise requires some specific assumptions about the pro- 
toplanetary disk (e.g., viscosity, scale height) which are 
poorly or not at all constrained empirically. Additionally, 
not knowing how empty the gap is in the disk mid-plane 
complicates any attempt to derive a useful companion 
mass estimate analytically and we leave it out of the 
scope of this letter. However, a direct search for com- 
panions could be attempted with high-contrast imaging 
instruments even though this might be challenging if the 
planet is stiU embedded (Wolf & D'Angelo 2005; Wolf 
et al. 2007). 

We strongly emphasize that there are alternative ex- 
planations for the observed disk morphology without in- 
voking an embedded planet. For instance, the disk rim 
could cast some shadow onto the disk surface. Puffed- 
up disk rims have been predicted at least for the inner 
edges of 'typical' protoplanetary disks even if the exact 
geometry and the question whether these rims can cause 
disk shadows is still debated (see, Dullemond & Mon- 
nier 2010, for a review). Concerning HD169142, the disk 
model by Honda ct al. (2012) docs not predict a shadow- 
ing disk rim. However, a combination of radially varying 
dust grain properties with dust settling might possibly 
be able to create a shadow-like feature without the need 
of a local deficiency in the surface density. Dust contin- 
uum observations with sufficient spatial resolution can 
shed light on this open question. 

4.3. The disk surface brightness 

The single power-law index we derive from fitting the 
surface brightness profile agrees with the results from 
Fukagawaet al. (2010), who found S oc 7--3-4±o.4 be^.^ygen 
120 - 200 AU in H band scattered light images. Grady 
et al. (2007) found a slightly shallower profile with S cx 
^-3.o±b.i between ~70 - 200 AU using HST/NICMOS 
at l.lyum. 



The errors in our profile fits justify the use of a broken 
power law. Broken power-law surface brightness pro- 
files are observed in resolved scattered light images of 
debris disks, where the inner disk regions show a shal- 
lower profile due to the presence of larger grains while 
the outer disk profile is steeper and dominated by blown- 
out smaller grains (e.g., Augereau ct al. 2001; Liu 2004). 
The disk around HD169142 is certainly still gas rich, but 
the observed brightness profile could be indicative of dif- 
ferent dust grain populations. However, also changes in 
the dust density distribution could lead to the observed 
effect. Until now, all studies focusing on the gas content 
or the dust contimmm of the HD169142 disk did not have 
enough spatial resolution to securely detect gaps or dis- 
continuities in the disk structure. Of particular interest 
would be to see whether the broken power law discovered 
here is accompanied with a radial change in the dust-to- 
gas ratio. 

Comparing the disk surface brightness derived by Fuk- 
agawa et al. (2010) with our values between 120 and 200 
AU suggests that roughly ~40% of the scattered light of 
HD169142 seems to be polarized in the H band. How- 
ever, as the errors bars in both studies are large, this 
value for the polarization fraction should be taken as a 
rough indication. At this level the polarization fraction 
appears to be higher than that for the HD100546 disk 
(^14%; Quanz et al. 2011) but comparable to the value 
for the transition disk around AB Aur (Perrin et al. 2009) 
for scattering angles around 90°. 

5. CONCLUSIONS 

Our PDI images of the protoplanetary disk around 
HD 169 142 reveal previously undetected features and em- 
phasize the power of this observing technique. With PDI 
on ground-based 8-m class telescopes the inner few tens 
of AU of the dusty surface layer of protoplanetary disks 
can be reliably resolved opening up a new era of disk 
imaging. Several PDI studies in the last few years have 
shown that sub-structures in disks (e.g., holes, gaps, spi- 
ral arms) appear to be rather the rule than the exception 
(e.g., Quanz et al. 2011; Hashimoto et al. 2011; Muto 
et al. 2012; Hashimoto et al. 2012; Mayama et al. 2012; 
Grady et al. 2013). 

In the case of HD169142, our images support previ- 
ous model results suggesting that the inner <20 AU are 
largely devoid of scattering dust particles. In addition, 
our images reveal an annular gap roughy between ^40 70 
AU, localized brightness asymmetries on opposite sides of 
this gap and a discontinuity in the radial surface bright- 
ness profile at ^120 AU. Our data alone do not allow 
us to unambiguously reveal the physical nature of the 
annular gap, i.c, whether it is a gap in surface density 
possible induced by forming planet, a shadowing efi'ect, 
or combination of radially changing dust properties and 
dust settling. However, our results make HD169142 an 
excellent target for follow-up studies aiming at a better 
understanding of disk physics and potentially ongoing 
planet formation. 

This research has made use of the SIMBAD database, 
operated at CDS, Strasbourg, France. We thank the 
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servations and M.R. Meyer, F. Meru and O. Panic for 
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